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A standard approach to investigate the non-perturbative QCD dynamics is through vacuum mod- 
els which emphasize the role played by specific gauge field fluctuations, such as instantons, monopoles 
or vortexes. The effective Hamiltonian describing the dynamics of the low-energy degrees of free- 
dom in such approaches is usually postulated phenomenologically, or obtained through uncontrolled 
approximations. In a recent paper, we have shown how lattice field theory simulations can be used 
to rigorously compute the effective Hamiltonian of arbitrary vacuum models by stochastically per- 
forming the path integral over all the vacuum field fluctuations which are not explicitly taken into 
account. In this work, we present the first illustrative application of such an approach to a gauge 
theory and we use it to compute the instanton size distribution in SU(2) gluon-dynamics in a fully 
model independent and parameter-free way. 



I. INTRODUCTION 



o 

Contemporary lattice gauge theory (LGT) simulations allow to compute from first principles a large class of hadronic 
matrix elements in QCD, in some cases even within a few percent accuracy. On the other hand, such simulations do 
not provide much detailed information about the structure of the gluonic fluctuations which drive the QCD dynamics 
in the strongly coupled regime. For example, despite several decades of investigations, the dynamical processes 
underlying chiral symmetry breaking and color confinement are still a matter of debate. 

The problem of identifying the dynamical origin of such non-perturbative phenomena has been extensively addressed 
O j in the context of phenomenological models which emphasize the role played by specific vacuum field fluctuations, such 
as e.g. instantons [I], monopololes [2] and center vortexes [3]. The configuration space of these models is defined by the 
^ H collective coordinates of the selected low-energy vacuum fields. On the other hand, the statistical distribution of such 
collective coordinates (or, equivalent ly, their effective Hamiltonian) is usually obtained through approximations upon 
which one does not have full theoretical control, e.g. by completely neglecting the contribution of the fluctuations 
around the chosen vacuum fields |4J or by estimating the role of such fluctuations through variational methods |5J. 

In principle, LGT simulations can be used to test the predictions of the phenomenological vacuum models, for 
example by looking for some specific signatures of the dynamics generated by instantons [U [7], monopoles [8] or 
vortexes |9J. On the other hand, the model dependence associated to the effective Hamiltonian for the vacuum field 
T—H degrees of freedom makes it is difficult to draw definitive conclusions about the validity of a given model. Indeed, 
a moderate disagreement with the experimental data or with the results of lattice QCD simulations may be due to 
X/-} either a wrong choice of the low-energy vacuum fields, or to the strong approximations involved in the definition of 
£^ their partition function. 

In order to tackle this problem, in a recent work we have developed a technique, which we shall refer to as Vacuum 
Manifold Projection (VMP), by which lattice simulations are used to rigorously compute the effective Hamiltonian of 
arbitrary vacuum models [10J, in a model independent way. This is done by non-perturbatively performing the path 
integral over all the vacuum field configurations which are not explicitly taken into account in the given low-energy 
vacuum model. For example, in an instanton model one performs the path integral over all the configurations which 
are orthogonal to the functional manifold spanned by multi-instanton configurations. Clearly, once the partition 
function has been evaluated from first principles, any failure of the model must be entirely due to the wrong choice 
of the effective low-energy degrees of freedom. 

In our first work, the VMP method was illustrated and tested by evaluating the instanton-antinstanton interaction 
in a simple quantum- mechanical toy model [TP] . Here, we present the first application to a gauge theory. In particular, 
we use the VMP method to compute the instanton size distribution in two-color Yang-Mills theory. 

The paper is organized as follows. In section IXT1 we review the VMP method for a generic choice of vacuum field 



degrees of freedom. In section III we focus on an effective theory based on instanton degrees of freedom and we use 
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the VMP method to compute the instanton size distribution in SU (2) gluon-dynamics. The main results, conclusions 



and perspectives are summarized in section IV 



II. THE VACUUM MANIFOLD PROJECTION METHOD 

Let us consider a gauge theory defined by the (Euclidean) path integral 

Z= ( VA^e~ s ^\ (1) 



where the S[A^\ formally includes the gauge- fixing and ghost terms, along with the fermionic determinant. In the 
following, we shall always assume that this path integral is defined in a Landau gauge. 

Let 7 = (71, . . . ,7/c) be a set of k collective coordinates which parametrize a manifold M. of vacuum field configu- 
rations A fI (x; 71, . . . , 7fc). For example, in instanton models, 71, . . . , 7jv are the positions, sizes and color orientations 
of all the pseudo-particles in the instanton ensemble. However, in general, we do not require the field configurations 
A^{x] 71, . . . , 7*;) to be solutions of the classical Yang-Mills equations of motion. 

For any given choice of the set of collective coordinates 7, a generic vacuum gauge field configuration A fI (x) 
contributing to the path integral can be decomposed as 

A^x) = i M (x;7i, . . . ,7fc) + B^x), (2) 

where B^(x) will be referred to as the "fluctuation field". Our goal is to use LGT to perform the path integral over 
such a field. More precisely, we want to compute the function ^(71, ... ,7/c) such that 

Z = J VA^e~ s[A ^ = Jd- fl ...d- fk (3) 

Eq. ([3| defines a statistical model in which 71, . . . , 7^ are the effective low-energy degrees of freedom and ^(71, ...,7/c) 
is the effective Hamiltonian. We shall see shortly that such function is defined as the logarithm of a gauge-fixed path 
integral. 

Since the representation (|3| of the path integral contains k additional integrals over dji, ^7^, we need to introduce 
k constraints. A natural choice is to impose a set of k orthogonality conditions: 

(5 M (x),^ 7i5/x (x,7)) ee Tr c |y > d 4 x 5 M (x) ^ 7i , M (x,7)| = 0, z = l,...,/c (4) 

9nA x ^) = ^r^te 7) ( 5 ) 

We observe that the functions {g 7 . (x, 7)}i=i,...,fc identify the k directions tangent to the manifold M. of background 
vacuum fields, in the point of curvilinear coordinates 7 = (71, . . . ,7^) — see Fig. [I] — . We consider only choices of 
the manifold M and of the point 7 such that the vectors ([5| define a system of coordinates. 

In the path integral formalism, the orthogonality conditions Q can be implemented by introducing a Fadeev-Popov 
representation of the unity. After some formal manipulation (see e.g. Ref. [5] [10] ) one arrives to an expression in the 
form of Eq. where the effective Hamiltonian ^(71, . . . , jk) is defined as 

«(7i,---,7fe) = -l°g{ J VBSfiM n*[(B M (x) )57( , M (x,7))] *[A(x, 7) + B(x)] e - S ^ x ^ +B ^ J, (6) 

where is a Jacobian factor and reads 

k 

^^KWl = [f[dji [ VU n 5[d,A% ]J5 - A M (x, 7 )^ 7llM (x,7))] ■ (7) 

J 1=1 J i 

U n (x) denotes a generic gauge transformation and A^(x) is result of gauge transforming the field A^(x) according to 
U Q (x). Note that, while the path integral Z is obviously gauge invariant, the definition of the effective Hamiltonian 
^(71, . . . , 7fc) relies on the choice of a Landau gauge. 
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Vacuum manifold 



Figure 1: Pictorical representation of the projection of the gauge field configuration A^{t) of the path integral Z — f VA^e~ s ^ A ^ 
onto a specific vacuum field manifold, spanned by two collective coordinates 71, and 72. A generic path is represented by a point 
in this three-dimensional space. The constraints given in Eq. Q imply that the the fluctuation field B^(x) is perpendicular to 
the plane tangent to the manifold in a given point. 



The orthogonality condition Q and the system of coordinates (5| can be used to device an algorithm to explicitly 
compute the effective Hamiltonian (J6|. We begin by observing that, on such a system of coordinates, a functional 
point A fI (x) which belongs to the vacuum field manifold M. is identified by the coordinates . . . , ^k) with 

*i[A M (a?)] = (i4 M (a?),(/ 74 , M (a?,7)). (8) 

Clearly, also field configurations which lie in a functional neighborhood of the manifold M. can be projected onto the 
same system of coordinates. In this case, for an arbitrary fixed choice of 7, the orthogonality conditions Q imply 
that there exits a set of collective coordinates 7 such that 

*i[A M (aO] = (A M (a;),(/ 7<iM (a:,7)) 

= (A^(x;j) + £ M (x),# 7i , /i (x,7)) 

= (A M (a?;7),0 7 . jM (aj,7)), (i = l,...,&)- (9) 

Such an equation allows to associate a set of collective coordinate 7 = (71, . . . , 7^) to any Landau gauge-fixed gluon 
field configuration A fI (x) which lies in the functional vicinity of the manifold A4. Based on this result, it is immediate 
to device an algorithm to evaluate the path integral 7/ (71, . . . , Jk) stochastically, from LGT simulations: 

1. An ensemble of N con f independent lattice configuations {[/^^(x), . . . , U^ Nconf \x)} is generated by standard 
LGT simulations. 

2. From such configurations, an ensemble of Landau gauge-fixed lattice configurations 
{U° . . . , E7* {Nconf \x)} is obtained, e.g. using the procedure illustrated in Ref.s [TT| [T2]. 

3. From each lattice configuration, the gluon field A 9 ^(x) at each lattice site is calculated from the gauge-fixed 
lattice link variables U^(x), for example using the discretized definition 

K{x) = ± (U°{x) - U°(x?) . (10) 

4. k non-linear equations for the 71, . . . ,7fe variables are obtained from Eq.s 

'(^(x),^ 1)M (x,7)) = (Af, fa-y),^ (x, 7)) 

(11) 

,{Al(x),gS k ^x^)) = (^(X7),<,>,7)). 
On the lattice, the inner products above are obviously represented as discretized sums, 

(Af l (x),g^Jx^)) = a A J2^c[Af 1 (n)g? /i (n;j)], (i = l,...,fe). (12) 
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where n runs over all the lattice sites and Trc refers to the trace over color labels. 



Notice that the quantities on the left- hand- side of Eq.s (11) are c- numbers which depend only on the lattice 
configuration and on the projection point 7. They correspond to the components of the lattice field A 9 ^(x) in 
the system of coordinates defined by the tangent vectors g^^x,^). 



On the other hand, the quantities on the right-hand-side of Eq. (11) are functions of the curvilinear coordinates 
7 which do not depend on the lattice configurations A 9 (x). These functions are completely specified, for any 
given fixed choice of the vacuum field manifold M and of the projection point 7. 



Hence, for any lattice configuration U 9 (x) the set of equations ( 11 ) can be numerically solved for 7 = (71, . . . , 7*.). 



5. The frequency histogram of 71, . . . , jk obtained by repeating this procedure for a large number N con f of statis- 
tically independent lattice gauge field configurations represents the probability V(ji, • . . , Jk) of a given set of 
curvilinear coordinates. Then the effective Hamiltonian for 71, ... ,7^ is obviously given by 

W( 7 i,..., 7fc ) = -logP(7i,...,7fc)- (13) 

Some comments on the VMP algorithm are in order. First of all, we emphasize that such a scheme relies on the 
assumption that the background field configurations A fI (x;j) are the relevant vacuum degrees of freedom. More 
precisely, we are requiring that the field configurations which contribute the most to the path integral lie in some 
functional neighborhood of the manifold A4 : so that they can be projected onto the system of coordinates defined 
in Eq. ([5]). Under such a condition, the definition of the effective Hamiltonian ^(71, . . . , jk) is unique and does not 
depend on any additional external parameters. 

Once the effective Hamiltonian has been determined, arbitrary vacuum-to- vacuum matrix elements can be approx- 
imatively computed by neglecting the contribution of the fluctuation field B^(x) to the operators. For example, if 
0{x) is a local operator which depends on the gluonic field A^ix), one has 

<0|d[A M (s)]|0> = \\ VA^O[A^(x)] e~ s ^ ~ ± fd 7l ...d 7k 0[i M (x; 7 )] e~ H ^-^ . (14) 

Since the effective Hamiltonian is evaluated from first principles, any violation of this identity implies that the 
fluctuation field plays an important role in the observable associated to this matrix elements or that there are 
important fluctuations in the path integral which are very far from the functional manifold .M, hence cannot be 
projected onto the system of coordinates In both cases, this would represent an unambiguous signature of the 
failure of the vacuum field model. 

We also stress that the VMP results do not depend on the choice of the gauge. Indeed, in Ref. [10J it was shown 



that once the system of Eq.s (11) is satisfied in one gauge, it holds also in any other gauge. 

Finally, we note that the VMP procedure is conceptually analog to a technique routinely adopted in classical 
statistical mechanics to evaluate the potential of mean- force G(ji : . . . , 7*.) as a function of a set of collective coordinates 
7ij ■ ■ • j7fc ( see e -g- H3] and references therein). In the canonical ensemble, the potential of mean- force (free energy) 
is defined as 



-/3G( 7 i,...,7fc) 



k 

JdrflSljt- fi(T)] 08 = l/k B T), (15) 



where T is the phase-space variable, H(T) is the Hamiltonian, and the functions fi(T) specify the definition of the 
macroscopic collective coordinate ji in terms of the microscopic phase- space variable T. The system's partition 
function reads 

Z = J d7i...d7 fc e -/^(7iv..,7 fc ) 5 ( 16 ) 

in complete analogy with Eq. ([3]). In order to evaluate G(7i, . . . , Jk) one generates and ensemble of statistically in- 
dependent equilibrium configurations {Ti, . . . , ^N conf }^ by means of Monte Carlo or molecular dynamics simulations, 
and evaluates 71, . . . ,7fc from 7^ = fi(T), for each of such configurations. In the limit of large number of equilib- 
rium configurations N con f the frequency histogram for the values 71 , . . . , jk obtained this way yields the equilibrium 
probability function for the collective coordinates, 7^(71, ... ,7fc). The negative of the logarithm of such a probability 
defines by construction the free energy 



-/3G(7i,..,,7fc) = logP(7i,...,7fc). 



(17) 
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III. CALCULATION OF THE INSTANTON SIZE DISTRIBUTION IN SU{2) GLUON- DYNAMICS. 

As a first illustrative application of the VMP method, in this section we specialize to the case in which the gauge 
theory is SU (2) gluon- dynamics and the vacuum manifold M is constructed by superimposing singular-gauge instanton 
and anti-instanton configurations, i.e. in the so-called instanton vacuum. 

Instanton vacuum models have been successfully used to investigate the dynamics of light hadrons [lj and the 
breaking of chiral symmetry [14] in QCD. It has been shown that this model is able to reproduce the spectrum of the 
observed lowest-lying hadrons [T5UTT] and of scalar and pseudo-scalar glueballs[18 - 21j, along with the electromagnetic 
structure of pions and nucleons [22-25J. In addition, the intanton liquid model provides an explanation of the AI = 1/2 
rule for non-leptonic hadron decays of kaons [26J and hyperons [27J, by promoting non-perturbative scalar diquark 
correlations [28J. While it is well known that semi-classically inspired vacuum models based only on singular-gauge 
instantons do not provide color confinement, a finite string tension was obtained in models employing cocktails of 
regular- and singular- gauge instantons and merons [4| [29] . 

From the theoretical point of view, a main limitation of the instanton approach to QCD resides in the well known 
"infrared catastrophe" of the semi-classical dilute instanton gas approximation: in the presence of quantum fluctu- 
ations, isolated instantons tend to swell. The instanton liquid model originates from the observation that such an 
infrared divergence can be cured if correlations between pseudo-particles are included [5j|35]. Unfortunately allowing 
for such interactions implies giving up a rigorous semi-classical theory of the QCD vacuum. In the so-called instan- 
ton liquid models the size distribution of the pseudo-particles is estimated using variational or phenomenological 
arguments. In such approaches, the typical instanton size in QCD is found to be of the order of 1/3 fm — see e.g. 

mm-- 

In order to extract information about the structure of instanton ensemble directly from LGT simulations, algorithms 
such as cooling [34] or eigenvalue filtering [7] have been proposed. Such methods provide techniques to filter out ( not 
integrate out) the high-frequency quantum content of lattice configurations. A problem with these methods is that 
their results critically depend on the choice of additional uncontrolled parameters, such as the number of cooling steps 
or of the number of retained low- lying eigen- modes of the Dirac operator. Clearly, such a dependence introduces some 
degree of arbitrariness in the results. For example, the instanton density vanishes if the limit of very large number of 
cooling steps. 

The goal of this section is to show that the VMP method can be used to rigorously evaluate the instanton size 
distribution directly from lattice simulations, without introducing any arbitrary parameter. The choice of focusing on 
SU (2) Yang-Mills theory was made in order to keep the analytical and numerical calculations as simple as possible. 



A. The Instanton Vacuum Manifold 



We begin our calculation by defining the manifold of vacuum field configurations. In two-color Yang-Mills theory, the 
classical field of an individual instanton or anti-instanton is specified by four collective coordinates: 7 = (z, p, #1, #2, #3) 
where z is the position of the pseudo-particle, p is its size, and 6i (i = 1, 2, 3) are three angles which specify a SU{2) 
matrix in color space according to: 



U = exp(iO k r k ). 



(18) 



A configuration of iV+ instanton and iV_ instantons A M (x;7i, 7at + +at_ ) can be constructed in the so-called sum 
ansatz, i.e. by superposing the classical fields of the pseudo-particles: 



iV_|_+iV_ 



(19) 



i=l 



where A 1 (x;ji) is the classical field of the i-th pseudo-particle: 



i*(x; 7i ) = t/V a C/t,y 



(x - z l ) v 



(X - Zi) 2 (X - Z*) 2 + p 2 ' 



(x - z l ) u 



(x - z*) 2 (x - z*) 2 + p 2 ' 



(for instantons) 
(for anti-instantons). 



(20) 
(21) 



77" and rj a are the so-called 't Hooft indexes. It is important to emphasize that the field obtained from the sum 
ansatz (19) is not in general a solution of Euclidean Yang-Mills equation of motion. 
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An effective theory for the two-color Yang-Mills theory can be obtained if the path-integral over all the gauge field 
configurations is replaced by the sum over all the configurations of a grand-canonical statistical ensemble of singular 
gauge instantons and antiinstantons: 

Z YM = £ ^^ c *(" + -"-> / n ^e-*(— -) (22) 

where 9 is the angle associated to strong CP violation and . . . , 7at) is the effective Hamiltonian, which represents 

the functional integral over the configurations which do not belong to the instanton vacuum manifold. In the analogy 
with classical statistical mechanics discussed in section [TTJ this term can be interpreted as the "potential of mean- force" 
between the pseudo-particles generated by all other quantum gauge field fluctuations in the path integral. Clearly, if 



the effective Hamiltonian W(7i, 7at) is calculated non-perturbativelly from first principles, then Eq.(22) provides 
an exact representation of the path-integral, which does not rely at all on the semi-classical approximation. 

If the instanton ensemble is not too dense, it is possible to perform a many-body expansion of the effective Hamil- 
tonian H(ji, ...,7at): 

H(7i,...,7at) ~ 5^ui(7i) + 5Z W2 ( 7i ' 7 ^ + Yl ^3(7i,7j,7fe) + ... (23) 

i i<j i<j<k 

The gauge invariance and translational invariance of the vacuum imply that the one-body term of the i— th pseudo- 
particle, 1/1(7*), is only a function of its size pi. Hence, the function n(pi) = exp[— 1/1(7*)] * s called the (single-) 
instanton size distribution. The divergence of such a term in the semi-classical dilute gas limit gives raise to the 
"infrared catastrophe". The terms 1^2 (7* > 7j)> ^3(7^ 7j> Ik): • • • describe multi-body interactions, and depend in general 
also on the relative positions and color orientations of the pseudo-particles. With such a definition, and for = the 
partition function reads 

Zym = E ^TaTT /II** II \ e-E*<i^(7* > 7i)+E*< i < fc t.3(7.,7i,7*)+». (24) 

N+,N- +' ~" i \ i=l J 

In the following, we present a computation of the single instanton size distribution n(p). The calculation of the 
many-body terms is conceptually analog and is referred to future work. 

B. Computing the instanton size distribution n(p). 

In order to calculate n(p), we have implemented the VMP method on the functional manifold spanned by individual 
instantons (and anti-instantons) defined in the singular gauge, which is a Landau gauge. The system of coordinates 
{g Jufl (x, 7)} was obtained by numerically differentiating the instanton field (20) and (21) with respect to each of the 
eight collective coordinates: 



g p ,v(x,p,z,{0i}) 



A^(x; p + Ap) - A^(x] p - Ap) 
2Ap 



9z»A x iPi z >{0i}) = — 2AzU » = 1,2,3,4) 

/ --snx\ Mx-6 l + A6)-A^x-6 l -A6) 
9eiA x >P> z >{0i}) = — 2A(9 , (z = 1,2,3), (25) 

where p, z 1 , z 2 , z 3 , z 4 , 61, 62, O3 define an arbitrary point on the manifold, whose choice will be specified below. 

In the VMP method, the distribution of collective coordinates is obtained by projecting lattice field theory config- 
urations, and solving the system of equations We used three different ensembles of SU(2) Landau gauge-fixed 
lattice QCD configurations, generated using the Wilson action [30] 1 : 

• an ensemble of 100 configurations generated on a 16 4 lattice at j3 = 2.3 (corresponding to a lattice spacing 
a ~ 0.17 fm) 



1 We thank Andre Sternbeck for providing us with such configurations. 
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VMP(a=0.17fm) 
VMP(a=0.12fm) 
VMP (a=0.09fm) 



-▼- Shuryak [95] 

Michael, Spencer [95] 



w 

WW 

m 



T. 



**w 



Instanton Size p(fm) 



Instanton Size p(fm) 



Figure 2: Left panel: Instanton size distribution n(p) obtained with the VMP technique, from three different sets of gauge- fixed 
configurations. The results have been obtained with p = 5a. Right panel: Results for pure gauge 577(2) obtained by C.Michael 
and S.Spencer from lattice simulations using the cooling algorithm (squares) |31| and by Shuryak from interacting instanton 
liquid model simulations (triangles) [32] . 



• an ensemble of 100 configurations generated on a 24 4 lattice at j3 = 2.4 (corresponding to a lattice spacing 



a ~ 0.12 fm) 

• an ensemble of 25 configurations generated on a 32 4 lattice at j3 = 2.5 (corresponding to a lattice spacing 
a ~ 0.09 fm). 

The gauge fixing of such configurations was performed using an over-relaxation algorithm — see [TTJ [12] — . 

For each lattice configuration we have used 10 different projection points 7 = (p, z 1 , z 2 , z 3 , z 4 , #1, #2? #3)- The color 
orientation angles §i and the space-time coordinate z v of the projection points were chosen randomly. On the other 
hand, the instanton size p of the projection point was held fixed to five lattice spacings, p = 5 a. This choice was 
made in order to ensure that p was always much larger than the lattice spacing a, yet much smaller than the size of 
the simulation box. 



Once the projection equations ( 11 ) were solved numerically for all lattice configurations and for all projection points, 
the probabilities for the collective coordinates p, z u , #i,02 and #3 were inferred from the corresponding frequency 
histograms. The distribution of the positions z, and color orientation angles 0i,#2>03 are trivial, as a consequence of 
the gauge invariance and translational invariance of the vacuum. On the other hand, the probability distribution n(p) 
for the instanton size p was found to display non-trivial structure, as it is expected as a consequence of dimensional 
transmutation. 

The instanton size distributions n(p) obtained from the three different ensembles of lattice configurations are shown 
in Fig. [2] We see that, within the statistical errors, the results are independent from the value of the lattice spacing a. 
The obtained distributions are peaked around p ~ 0.2 fm and qualitatively agree with the results of other methods 
(see right panel of Fig. |2|. On the other hand, we stress that the present results do not depend on any arbitrary 
parameter, and take into account the full quantum content of the lattice configurations. 

It is instructive to compare the small instanton size tail of the n(p) distribution obtained non-perturbatively through 
the VMP method with the analytical one-loop formula obtained by 't Hooft [33J, 

n ne-loop(p) OC p^ Nc ~ 5 . (26) 

In Fig. [3] we show that our results are in quantitative agreement with the leading-order perturbative prediction in 
the range 0.1 fm < p < 0.2 fm. The VMP results for very small-sized instantons (p < 0.1 fm) are affected by lattice 
discretization errors, while the perturbative calculation is not reliable for large instanton sizes, p > 1/Asu(2)- Such a 
comparison shows that the suppression of the large-sized instantons is a purely non-perturbative effect. 



C. Testing the accuracy of the VMP calculation of the instanton size 



In this section we present a study of the accuracy of our VMP calculation of the instanton size distribution based on 
an instanton model. We generated 1000 configurations of an ensemble of configurations constructed by superimposing 
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Figure 3: Logarithm of the instanton size distribution n(p), obtained with the VMP technique from two different sets of gauge- 
fixed configurations. The results are compared with the exact one- loop calculation obtained by 't Hooft [33]. The VMP results 
are in agreement with the one-loop calculation, in the small-sized instantons range 0.1 fm < p < 0.2 fm. The VMP results for 
very small-sized instantons (p < 0.1 fm) is affected by lattice discretization errors. 



the fields of 20 instantons and 20 anti-instantons in a box of volume V = (2.7 fm) 4 . The positions and color orientations 
of the pseudo-particles were randomly chosen, while the size of the pseudo-particles was sampled from a Gaussian 
distribution. 



^model 

(p) oc exp 



(P-Po) 2 

2<7 2 



(0 fm < p < 1 fm), 



(27) 



with po = 0.3 fm, and a = 0.13 fm. The field configurations were discretized on a lattice of spacing a = 0.17 fm. The 
VMP procedure outlined above was applied to such an ensemble of instanton-model lattice configurations in order to 
compute instanton size distribution n(p). In Fig. [4] we compare the distribution n(p) calculated through the VMP 
method and the exact n mo dei(p) which was used to generate the instantons ensemble. We see that the VMP approach 
allows to very accurately reconstruct the correct instanton size distribution for all values of the instanton size larger 
than about a lattice spacing. 



IV. CONCLUSIONS 



In this work we have presented the first application of the recently developed VMP method to a gauge theory. 
We have used such a method to perform a non-perturbative calculation of the instanton size distribution in SU (2) 
gluon-dynamics. Our results are in good agreement with other lattice calculations for all p, and with the one- loop 
perturbative estimate, in the small instanton size regime. The VMP calculations of the many-body terms of the 
effective Hamiltonian of the instanton vacuum models is conceptually analog to the one presented here. For example, 
in order to compute the instanton-instanton two-body interaction i^Ti, Ik) one would need to project onto the 
manifold spanned by the collective coordinates of two pseudo-particles. Such a calculation was performed in [10] in 
the case of a quantum mechanical toy model and its extension to gauge theories does not raise conceptual difficulties. 

It should be emphasized that, in the present calculation, the fluctuations around the instanton field are not just 
filtered away, as in the cooling or eigenvalue filtering algorithms. Instead, they are systematically integrated out. As 
a result, the calculation is entirely self-consistent and does not depend on any arbitrary external parameter. Recently, 
Perez and co-workers have proposed a so-called Adjoint Filtering Method based on Dirac quasi-zero modes in the 
adjoint representation [36J, which does not rely on additional parameters. This method has been used to study 
the dynamics of smooth single-instanton configurations "heated" by Monte Carlo updating steps. The present VMP 
method discussed in this paper could be used in principle to perform the same task. On the other hand, the VMP 
approach can be applied also to vacuum field configurations which are not solutions of the equation of motion. 




Figure 4: Comparison between the distribution used to generate an instanton ensemble with density n ~ 0.75(/ra 4 ) and the 
distribution obtained using the VMP approach from the corresponding lattice configurations. 
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